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ABSTRACT 

Thick thermal barrier coating systems in a diesel engine experience severe thermal 
low cycle fatigue (LCF) and high cycle fatigue (HCF) during engine operation. In the 
present study, the mechanisms of fatigue crack initiation and propagation, as well as of 
coating failure, under thermal loads which simulate engine conditions, are investigated 
using a high power CO 2 laser. In general, surface vertical cracks initiate early and grow 
continuously under LCF and HCF cyclic stresses. It is found that in the absence of 
interfacial oxidation, the failure associated with LCF is closely related to coating sintering 
and creep at high temperatures, which induce tensile stresses in the coating after cooling. 
Experiments show that the HCF cycles are very damaging to the coating systems. The 
combined LCF and HCF tests produced more severe coating surface cracking, 
microspallation and accelerated crack growth, as compared to the pure LCF test. It is 
suggested that the HCF component cannot only accelerate the surface crack initiation, but 
also interact with the LCF by contributing to the crack growth at high temperatures. The 
increased LCF stress intensity at the crack tip due to the HCF component enhances the 
subsequent LCF crack growth. Conversely, since a faster HCF crack growth rate will be 
expected with lower effective compressive stresses in the coating, the LCF cycles also 
facilitate the HCF crack growth at high temperatures by stress relaxation process. A surface 
wedging model has been proposed to account for the HCF crack growth in the coating 
system. This mechanism predicts that HCF damage effect increases with increasing 
temperature swing, the thermal expansion coefficient and the elastic modulus of the 
ceramic coating, as well as the HCF interacting depth. A good agreement has been found 
between the analysis and experimental evidence. 


t National Research Council — NASA Research Associate at Lewis Research Center. 
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INTRODUCTION 


Ceramic thermal barrier coatings have attracted increasing attention in heat engines 
because of their ability to provide thermal insulation to engine components. The advantages 
of using the ceramic coatings include a potential increase in engine operating temperature 
with elimination of the water cooling system and a longer service life in the harsh in- 
cylinder environment. Zr02-based ceramics are the most important coating materials for 
such applications because of their low thermal conductivity, relatively high thermal 
expansivity and excellent mechanical properties. A typical thermal barrier coating system 
consists of a top layer Zr02-8%Y203 coating and an intermediate superalloy-type bond 
coat and the alloy substrate. The application of advanced thick thermal barrier coatings 
(TTBCs) for diesel engine components such as piston crowns and cylinder heads is 
promising for increasing engine fuel efficiency, performance and reliability [1, l 

However, durability of thick thermal barrier coatings under severe temperature 
cycling conditions encountered in a diesel engine remains a major problem. In a diesel 
engine, two types of thermal fatigue transients exist [1, 3 ’ 4] . The first transient type, which 
is associated with the start/stop and no-load/full-load engine cycle, generates thermal low 
cycle fatigue (LCF) in the coating system. The second transient type, which is associated 
with the in-cylinder combustion process, generates a thermal high cycle fatigue (HCF). It 
occurs at a frequency on the order of 10 Hz (i.e., 1000-2600 RPM). The HCF transient can 
generate a temperature fluctuation of more than 200°C that will superimpose onto the 
steady-state engine temperature at the coating surface 3 ’ 5 l Therefore, the failure 
mechanisms of thick thermal barrier coatings are expected to be quite different from those 
of thin TBCs under these temperature transients. The coating failure is related not only to 
thermal expansion mismatch and oxidation of the bond coats and substrates ^ 6 ’ 1 \ but 
also to the steep thermal stress gradients induced from the temperature distributions during 
the thermal transients in the coating systems [1, 2 ’ 7-101 . 

The development of advanced thick thermal barrier coatings requires a thorough 
understanding of thermal fatigue behavior. Although it has been reported ^ 111 that stresses 
generated by a thermal transient can initiate surface and interface cracks in a coating 
system, the mechanisms of the crack propagation and of coating failure under the complex 
LCF and HCF conditions are still not understood. Particularly, the understanding of surface 
vertical crack propagation in thick thermal barrier coatings under thermal cyclic loading is 
of great importance. Experimental evidence has shown all coating failure under severe 
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thermal cycling conditions, produced either by a high heat flux burner rig or a high power 
laser, is more or less associated with surface vertical cracks . These vertical surface 
cracks and sometimes through-thickness-cracks can facilitate the interfacial crack 
formation, eventually resulting in the coating delamination and spallation. In addition, the 
interaction between LCF and HCF cycles, and the impact of relative amplitude of the LCF 
and HCF transients on coating fatigue life are among the most important aspects in 
understanding the thermal fatigue behavior of the coating systems. In this paper, thermal 
fatigue behavior of an yttria partially stabilized zirconia coating system under simulated 
LCF and HCF engine conditions is investigated. The effects of LCF and HCF parameters 
on surface fatigue crack initiation and propagation in the coating are also discussed. 

EXPERIMENTAL MATERIALS AND METHODS 
Materials and Specimen Preparation 

Zr02-8 wt % Y2O3 ceramic coating and Fe-25Cr-5Al-0.5Y bond coat were 
plasma-sprayed onto 4140 and 1020 steel substrates using an ABB ASEA IFB2000 6-axis 
industrial robot. The plasma spray conditions used for both the ceramic coating and bond 
coat are listed in Table 1. The sample substrate configurations were rectangular bar, as well 
as angle iron which provided a comer shape for the coating. The specimen dimensions are 
illustrated in Figure 1. The thickness of the ceramic coating was about 1.5- 1.6 mm. The 
bond coat thicknesses were 0.28 mm and 0.5 mm for the angle iron specimens and the 
rectangular flat specimens, respectively. 

Table 1 Plasma spray parameters for ZrQ2-8wt%Y2p3 top coat and FeCrAlY bond coat 


Coatings 

materials 

Torch 

power 

KW 


Carrier 
gas flow 

Standard 

liter/min. 

Spray 

distance 

mm 

Feed 

rate 

g/min. 

Torch 

translation 

rate 

mm/s 

Air 

cooling 

condition 

Psi 

Substrate 

temperature 

°C 

FeCrAlY 
PRAX- 
AIR 
FE213 
44-74 pm 

35 

(9mB 
plasma 
torch, GH 
nozzle) 

Ar: 56.6 
N 2 : 9.4 

Ar: 8.3 

127 

68 

1300 

50 

250 

Zr0 2 - 
8 %Y 2 03 
ZIRCOA 
9507/46 
44-74 pm 

40 

(9mB 
plasma 
torch, GH 
nozzle) 

Ar: 14.2 
N 2 : 7.1 

Ar: 3.2 

101.6 

20 

1000 

50 

250 
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Low Cycle and High Cycle Fatigue Tests 

Low cycle and high cycle fatigue tests under simulated engine temperature and 
stress conditions were conducted using a high power 1.5 KW CO 2 laser (EVERLASE, 
Coherent General Inc., Massachusetts). This test rig was controlled by a PC programmed 
to simulate different LCF and HCF temperature cycles. In this study, the HCF combustion 
cycles were simulated using the pulsed laser mode. The laser pulse period and pulse width 
were set at 92 and 9 milliseconds (ms) respectively, with effective square wave equivalent 
pulse heating time about 6 ms. The total beam power in the pulsed mode was set to 
approximately 180W. The laser pulse input waveform, measured by an oscilloscope (THS 
720 Tekscope with frequency 100 MHz and data acquisition rate 500 Meg samples/sec., 
Tektronix, Oregon), is shown in Figure 2. 


Laser pulse waveform 



Time, ms 


Fig. 2 Laser pulse waveform recorded from the laser pulse signal by THS 720 
Tekscope. 


The laser power density for an idealized spherical Gaussian beam is related to laser 
total power P and beam radius w by the following relation I13, 14] 
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where 7 0 is laser power density at the center, r is the distance from the center. The beam 
radius w has been defined as the distance at which the laser power density has dropped to 
1/e 2 of its value at the center. In this study, in order to produce a lower power density 
suitable for simulating diesel engine conditions, and also to cover a larger test specimen 
area, a Plano Concave ZeSe lens with focal length -330 mm was used to expand the laser 
beam. With the specimen being placed at a distance 460 mm from the magnifying lens, the 
beam radius w was increased from 7 mm to about 16 mm, as determined from laser bum 
patterns. Laser power density distributions under the test conditions are shown in Figure 3. 


Expanded beam, pulse mode Raw beam, pulse mode 

Expanded beam, CW mode 



Fig. 3 Laser power density distributions estimated from the measured laser waveform 
and total power output. Minor beam non-uniformity observed is neglected. 


During the thermal fatigue testing, specimen surface temperatures were measured 
by two 8 micron infrared Pyrometers (Model MX-M803 Maxline Infrared Thermometer 
Measurement and Control System, Ircon, Inc., Illinois), aimed at the beam center (giving 
the peak temperature) and 7 mm away from the center, as shown in Figure 1. The backside 
metal temperature was determined by an R-type thermocouple. For the combined LCF and 
HCF tests, the pulsed laser mode was used to generate the heating and cooling cycles, and 
the total power output was 180W. Two sets of experiments were conducted for angle iron 
specimens, with heating/cooling cycle times set at 30/5 and 5/3 minutes respectively. 
Because the high energy laser pulse was used, an HCF component was inherently 
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superimposed on the LCF cycles. These experiments were designed to provide information 
on LCF and HCF interactions, and the effect of relative LCF and HCF cycle numbers on 
ceramic coating failure mechanisms. Backside air cooling was used to maintain the desired 
temperature gradient. The backside metal temperature was fixed at about 250°C, by simply 
adjusting the cooling air flow. Steady state heating was usually reached in two to three 
minutes. The peak specimen surface temperature (steady-state average temperature at the 
beam center location) thus measured was about 850°C. The total HCF cycle numbers were 
fixed at 10 x 10 6 cycles for the angle iron specimens, corresponding to a total heating time 
of about 256 hours. In order to study the effect of surface temperature on fatigue behavior, 
another angle iron test was conducted of using a backside temperature fixed at 350°C, with 
a corresponding surface center temperature about 950°C. A pure LCF test was also 
conducted using the continuous wave (CW) laser, with a same total power 180W and a 5 
minute heating and 3 minute cooling cycle, to study coating fatigue behavior in the absence 
of an HCF component. A similar set of pure LCF and combined LCF and HCF tests were 
also carried out for the rectangular flat specimens. With a fixed back temperature of 250°C, 
the 180W pulsed laser beam generated a surface center temperature of approximately 
920°C. The tests were used to provide information on crack distributions and coating 
fatigue behavior of flat specimens. The specimen and experimental conditions for LCF and 
HCF tests are summarized in Table 2. 


Table 2. Summary of specimen and experimental conditions 


No 

Material 

Test type 

Surface 

tempera- 

ture 

°C 

Backside 
metal 
tempera- 
ture °C 

Heating/ 

cooling 

time, 

min. 

Total 

heating 

time 

hrs. 

Total 

HCF 

cycles 

xlO 6 

Total 

LCF 

cycles 

1 

Angle iron TBC 
tc= 1.6mm 
tb=0.28mm 

LCF 

CW 180W 

850 

250 

5/3 

256 

- 

3067 

2 

Angle iron TBC 
tc= 1.6mm 
tb=0.28mm 

LCF& 

HCF 

Pulse 180W 

850 

250 

30/5 

256 

10 

510 

3 

Angle iron TBC 
tc= 1.6mm 
tb=0.28mm 

LCF& 

HCF 

Pulse 180W 

850 

250 

5/3 

256 

10 

3067 

4 

Angle iron TBC 
tc= 1.6mm 
tb=0.28mm 

LCF& 

HCF 

Pulse 180W 

950 

350 

30/6 

256 

10 

510 

5 

Flat TBC 
tc= 1.5mm 
tb=0.5mm 

LCF& 

HCF 

Pulse 180W 

920 

250 

30/5 

153 

6 

307 

6 

Flat TBC 
tc= 1.5mm 
tb=0.5mm 

LCF 

CW 180W 

920 

250 

30/5 

153 

- 

307 
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Since the pyrometer has a slower response time (> 25 ms) compared to the actual 
laser pulse width (6 ms), the temperature swing generated by the pulsed laser on the 
ceramic surface could not be recorded. Therefore, one dimensional finite difference 
analysis has been used to model the thermal HCF temperature profile, providing the 
important thermal parameters such as the temperature fluctuation AT and interaction depth 
on the ceramic surface under the given test conditions. 

Microscopic Examinations 

The tested coating surfaces and cross-sections were examined under both optical 
and electron scanning microscopes to obtain information on crack density and distribution, 
as well as crack surface morphology. To prevent damage by specimen cross-section 
preparation, a pressurized epoxy infiltration method for specimen mounting was used. By 
this technique, epoxy was first poured over the specimens and their holding cups in a 
vacuum chamber. After the epoxy degassing in vacuum, the specimens were moved into a 
pressurized chamber (up to 1200 Psi) for 24 hours, as the epoxy cured. Therefore, the 
epoxy filled the cracks in the specimen, and the original crack characteristics generated in 
thermal fatigue tests were preserved. 

EXPERIMENTAL RESULTS 

Temperature Cycles Induced by Laser Beam Heating 

Figure 4 shows typical temperature cycles of laser thermal fatigue tests. The steady 
states were reached during the first few minutes of the cycling. It may be noticed that under 
the combined LCF and HCF conditions, even though the pyrometer could not accurately 
read the temperature fluctuations of the HCF component because of its slow response time, 
large variations in recorded temperatures were still observed during laser heating. In 
contrast, the continuous wave laser test simulating the pure LCF condition showed very 
little temperature fluctuation. This suggests that regardless of the similar steady state 
average temperature profiles produced by the pulsed laser beam and the CW laser beam, 
the pulsed laser beam heating induced a severe surface temperature swing which was 
superimposed onto the steady state temperature. 

Because of an expanded near-Gaussian laser beam used, temperature distributions 
are expected to vary across the beam diameter. This was confirmed by experiments, as 
shown in Figure 4. The average temperature reading from the pyrometer aimed at a point 
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7 mm away from the center is 250°C lower than that from the pyrometer aimed at the 
center for the angle iron specimens. Even higher temperature differences were observed for 
the flat specimens. This Gaussian beam profile, in principle, can provide additional 
information on coating failure mechanisms with heat flux distributions, establishing a 
relationship between the coating damage and the test parameters, such as the average 
surface temperature and temperature swing from a set of experiments. 

Temperature and Thermal Stress Distributions 

Figure 5 shows the calculated temperature distributions (with a simplified one- 
dimensional configuration) across the thermal barrier coating system on an angle iron 
during the steady state heating under various heat fluxes. Because of the constraints 
imposed by the angle iron structure, specimen bending was not likely to occur. Therefore, 
the in-plane stress distributions in the system at the steady state during the first heat up 
could be calculated from the mechanical equilibrium and strain compatibility conditions. 
The results are shown in Figure 6. The material properties used in the calculations are listed 
in Table 3. It should be noted that the overall stress is the summation of the thermal stress 
and residual stress in the system. As will be discussed later, for longer heating times, 
ceramic sintering and creep will become significant, thus modifying the stress states in the 
coating system. 

When pulsed laser heating is used, a severe thermal transient will be induced even 
in the absence of LCF cycling. This temperature fluctuation and history under the HCF 
conditions were modeled by the one dimensional finite difference approach. In order to 
verify the validity of this model under the present laser beam conditions, the one 
dimensional finite difference analysis method was compared with analytical solutions for 
both a uniform, constant irradiance model and a Gaussian beam model in calculating the 
surface temperature swing ^ 12 l The temperature swing predicted by all three approaches 
was essentially the same, implying that the Gaussian beam is sufficiently widespread to 
allow the use of the one-dimensional assumption. The modeled results indicate that the 
HCF transient occurs only at the surface layer of the ceramic coating. This layer may be 
defined as the HCF interaction depth at which appreciable temperature fluctuation (20°C or 
above) will occur. This temperature swing generated by the pulsed laser increases with 
increasing the laser peak power density and the laser pulse width (laser pulse heating time), 
as shown in Figure 7. However, the HCF interaction layer depth, which is independent of 
laser power density, increases with increasing laser pulse width, as illustrated in Figure 8. 
Under the HCF condition of 6 ms heating, the interaction depth is about 0.15 mm, as 
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calculated by the finite difference method. The HCF component, therefore, is generated 
only on the very surface of the ceramic coating. However, the effect of HCF on thermal 
fatigue is more complex and will extend far beyond this characteristic depth, as will be 
discussed later. 

The temperature profiles generated by the pulsed laser under peak heat fluxes 3.38 
and 4.95 MW/m 2 are illustrated in Figure 9. The HCF stress distributions with coating 
depth and variations with time are shown in Figure 10. It can be seen that this temperature 
fluctuation induces high-frequency cyclic stresses on the coating surface, with the predicted 
HCF stress ranging from around 60 MPa at 3.38 MW/m 2 to 100 MPa at 4.95 MW/m 2 . 
The dashed lines in Figure 10 represent the ceramic surface stress values at the average 
steady state surface temperatures under the corresponding average heat fluxes 0.220 and 
0.323 MW/m 2 , respectively. 


Table 3 Physical and mechanical properties of the thermal barrier coating system 


used in calculations 


Material Properties 

Plasma sprayed Zr0 2 - 
8%Y 2 0 3 

Plasma sprayed 
FeCrAlY 

Steel substrate 

Thermal conductivity 
k, W/m-K 

0.9 

11.0 

46.7 

Thermal expansion 
coefficient 
a , m/m °K 

10.8 X 10' 6 

12.4 XI O' 6 m/m°C 

14.2 XI O' 6 m/m°C 

Density 
p , kg/m 3 

5236 

- 

7850 

Heat capacity 
c, J/kg-K 

582 

- 

456.4 

Young's modulus 
E, GPa 

27.6 

137.9 

207.0 

Poisson's ratio, V 

0.25 

0.27 

0.25 
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Temperature, °C Temperature, 


(a) LCF (CW Laser) Thermal Fatigue Test 


(b) Combined LCF and HCF Thermal Fatigue Test 



Time, minutes 



(c) Combined LCF and HCF Thermal Fatigue Test 
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800.0 
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(d) Combined LCF and HCF Thermal Fatigue Test 
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JO) 


(e) Combined LCF and HCF Thermal Fatigue Test 



Fig. 4 Laser heating and cooling profiles simulating the engine operating conditions 
(temperature has been corrected with angle iron configurations), (a) Angle iron 
with LCF only; (b), (c) and (d) Angle iron with combined LCF and HCF; (e) Flat 
specimen with combined LCF and HCF. 
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Distance from the surface, mm 

Fig. 5 Temperature distributions in a thermal barrier coating system on an angle iron at 
steady state heating for various heat fluxes. 


-a - - 0 th , 0.10 MW/m 2 — • o re , ceramic a o t0tal , 0.10 MW/m 2 

□ - - 0 th , 0.20 MW/m 2 - - » c f e , bond coat □ G total , 0.20 MW/m 2 

o - - 0 th , 0.32 MW/m 2 - -a- - - substrate o a total , 0.32 MW/m 2 
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50.0 

£ 00 

§ 

-50.0 
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£ -100.0 
-150.0 
- 200.0 
-250.0 

0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 

Distance from the surface, mm 

Fig. 6 Stress distributions in a thermal barrier coating system on an angle iron at steady 
state heating for various heat fluxes. 
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Fig. 7 


Fig. 8 



Maximum temperature fluctuation, calculated by finite difference analysis, as a 
function of laser pulse width. 



Laser pulse width, ms 


The relationship between laser interaction depth and laser pulse width. The 
interaction depth is independent of laser power density. 
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Temperature, °C Temperature, 



Distance from the surface, mm 

(a) 



Distance from the surface, mm 


(b) 

Fig. 9 Predicted temperature profiles generated by pulsed laser heating (pulse width 6 
ms). A higher heat flux produces a higher temperature swing and average surface 
temperature (dashed line represents the average temperature). 
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Stress, MPa Stress, MPa 




Stress distributions in thermal barrier coating during thermal HCF test 
50.00 

0.00 

-50.00 

- 100.00 

-150.00 

- 200.00 

-250.00 

0.00 0.20 0.40 0.60 0.80 1.00 1.20 1.40 1.60 
Distance from the surface, mm 

(a) 

Stress Temperature 



Thermal HCF stress history at the ceramic surface 



Time, seconds 
(b) 

Fig. 10 Predicted thermal stresses induced by pulsed laser heating. Besides a constant 
stress gradient generated by the steady state heating, high frequency HCF cyclic 
stresses are present near the ceramic coating surface. Peak power density 3.38 
MW/m 2 . 
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Stress, MPa Stress, MPa 


Stress distributions in thermal barrier coating during thermal HCF test 



Distance from the surface, mm 

(c) 


Stress Temperature 


Fig. 10 


Thermal HCF stress history at the ceramic surface 
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Predicted thermal stresses induced by pulsed laser heating. Besides a constant 
stress gradient generated by the steady state heating, high frequency HCF cyclic 
stresses are present near the ceramic coating surface (continued). Peak power 
density 4.95 MW/m 2 . 
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LCF and HCF Damage on Thermal Barrier Coatings 

The surface cracking was observed for all specimens tested under LCF and/or HCF 
conditions (total heating time up to 256 hours). Compared to the pure LCF tested 
specimen, the combined LCF and HCF tests produced much higher crack densities, with 
more complex crack networks on the ceramic surfaces. Examination of surface cracks on 
the flat specimens shows that the crack density decreases with decreasing laser power 
density. 

The crack patterns on the angle iron and flat specimen surfaces are schematically 
illustrated in Figure 11. At the angle iron comers, nearly parallel cracks which ran across 
the comers were formed by the laser thermal fatigue tests. In contrast, equiaxial crack 
networks (mud flat cracks) were generated by the laser beam at the flat specimen surfaces. 
However, at the edges of the flat specimens, parallel cracks similar to those found on the 
angle iron comers were observed with crack direction perpendicular to the edges. 
Compared to pure LCF tests, the combined LCF and HCF initiated more secondary cracks, 
and micro- spallation at the cracked surfaces. The optical micrographs of the cracked 
surfaces are shown in Figure 12. The results suggest that much higher surface stresses 
were induced at the ceramic surface by the pulsed laser HCF component. 




(b) Rectangular flat specimen 

Fig. 11 Schematic diagram showing the crack patterns on coating surfaces after laser 
testing. 
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Figure 13 shows SEM micrographs of the tested coatings on angle iron specimens. 
It can be noticed that the pure LCF tested specimen shown in Figure 13 (a) has the most 
intact coating surface, and the thermal fatigue cracks are relatively regular with well 
matched crack faces. However, the combined LCF and HCF tests produced more severe 
coating surface damage. Besides the major thermal fatigue cracks, surface coating micro- 
spallation, crack branching and loose particles intruding into the cracks are often observed. 
For all combined LCF and HCF tested specimens, the specimen with the 30 minute 
heating/5 minute cooling cycles at a lower temperature (850°C) showed the least surface 
damage. In contrast, the most surface damage was found for the specimen with the 30 
minute heating/5 minute cooling cycles at the higher temperature (950°C). In the latter 
specimen, cracks were branched into multiple crack networks and accompanied with more 
coating spallation, and the major crack density and the crack width were also significantly 
higher compared with the lower temperature tested specimens. 

DISCUSSION 

Ceramic Coating Sintering and Creep at High Temperatures 

During thermal fatigue testing, ceramic sintering and creep will occur under the 
given temperature and stress conditions. Due to the porous and microcracked nature of 
plasma-sprayed ceramic coatings, the primary creep stage is often observed for these 
coatings, with the strain rate continuously decreasing with time ^ 15, 16 \ This creep behavior 
is probably related to stress-enhanced ceramic sintering phenomenon, the splat relative 
sliding, and the stress redistribution around the splats and microcracks. The stress- 

dependent deformation can result in coating shrinkage and thus stress relaxation at 
temperature under the compressive thermoelastic stresses. The strain rate £ p can be 

generally written as 


e p = A- exp- • (^thf S (2) 

where A , n and s are constants, Q is activation energy, R is gas constant, o th is the in- 
plane compressive thermal stress in the coating, and t is time. The time exponent s is 
reported to be 0.82 under low stresses (<80 MPa), and to be 0.67 under high stresses (up 
to 655 MPa) [15, 161 . The creep strain £ p l in the ceramic coating can be expressed as 
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£p = je p (cr Ih ,T,t)dt = jAexp 
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where e p l and e p l ~ l are creep strains at time t h and the previous time step t t _ x , 
respectively, is the initial thermal compressive stress in the coating, E c and v c are the 

elastic modulus and Poisson's ratio of the ceramic coating. The stress relaxation effect on 
the total creep strain is considered by the e p l ~ l term in Equation (3). Using the literature 

reported data A, n, s and Q for the plasma-sprayed ceramic coating [15, 16] , the creep 
strains as a function of time can be estimated for a heat flux 0.323 MW/m 2 , as illustrated in 
Figure 14 (a) and (c). The in-plane stress distribution profiles in the coating, as shown in 
Figure 14 (b) and (d), indicate that significant stress relaxation will occur, especially at the 
top half of the coating, because of higher thermal stresses and temperatures at these 
locations. In addition, the creep strain and thus stress relaxation increase with decreasing 
the time exponent 5 . The coating creep and stress relaxation are strongly dependent upon 
the stress exponent, n, and the activation energy, Q. As illustrated in Figures 14 and 15, 
with a higher n value and a slightly lower activation energy, more significant stress 
relaxation will occur in the coating system. 


The laser heat flux has a significant effect on coating creep and stress relaxation. As 
shown in Figure 16, a lower laser heat flux (0.20 MW/m 2 ) will establish a lower surface 
temperature and a less steep temperature gradient across the coating, therefore, a lower 
thermal stress will be expected in the coating. As a consequence, total creep strain and 
stress relaxation will be much less as compared with those in the high heat flux case. 
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(a) 



(b) 


Fig. 12 Optical micrographs showing the cracked coating surfaces after laser thermal 
fatigue testing, (a) and (b) The coating surface with pure LCF test; 
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(C) 



(d) 


Fig. 12 Optical micrographs showing the cracked coating surfaces after laser thermal 
fatigue testing (continued), (c) and (d) The coating surface with LCF+HCF 
test (arrows show regions with imminent spalling); 


NASA TP-3676 


21 




(e) 

Fig. 12 Optical micrographs showing the cracked coating surfaces after laser thermal 
fatigue testing (continued), (e) The coating edge with LCF+HCF test. 
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(a) 


Fig. 13 SEM micrographs showing the coating surface morphologies after laser 
thermal LCF and HCF testing for angle iron specimens, (a) LCF tested, 
5 min. heating/3 min. cooling cycle, center temperature 850°C. 
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(b) 

Fig. 13 SEM micrographs showing the coating surface morphologies after laser thermal 
LCF and HCF testing for angle iron specimens (continued), (b) LCF+HCF tested, 
30 min. heating/5 min. cooling cycle, center temperature 850°C. 
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(c) 

Fig. 13 SEM micrographs showing the coating surface morphologies after laser thermal 
LCF and HCF testing for angle iron specimens (continued), (c) LCF+HCF 
tested, 5 min. heating/3 min. cooling cycle, center temperature 850°C. 
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(d) 

Fig. 13 SEM micrographs showing the coating surface morphologies after laser thermal 
LCF and HCF testing for angle iron specimens (continued), (d) LCF+HCF 
tested, 30 min. heating/6 min. cooling cycle, center temperature 950°C. 
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Fig. 14 The creep strains and stress relaxation in the ceramic coating as a function of 
time. The results are estimated from available literature data for the case of heat 
flux 0.32 MW/m 2 . The total strains and stress relaxation at different layer depths 
in the ceramic coating increase with the time exponent s . (a) and (b) 5 = 0.82. 


NASA TP-3676 


27 




Creep strain, % 



Time, hours 

(c) 


Time, hours 



Distance from the surface, mm 

(d) 

Fig. 14 The creep strains and stress relaxation in the ceramic coating as a function of 
time. The results are estimated from available literature data for the case of heat 
flux 0.32 MW/m 2 . The total strains and stress relaxation at different layer depths 
in the ceramic coating increase with the time exponent s . (c) and (d) 5 = 0.67 . 
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(b) 

Fig. 15 The creep strains and stress relaxation in the ceramic coating as a function of 
time. Compared with Fig. 14, the total creep strains and stress relaxation in the 
ceramic coating are increased with a higher stress exponent n and a lower 
activation energy Q. (a) and (b) n = 0.8. 
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Fig. 15 The creep strains and stress relaxation in the ceramic coating as a function of 
time. Compared with Fig. 11, the total creep strains and stress relaxation in the 
ceramic coating are increased with a higher stress exponent n and a lower 
activation energy Q. (c) and (d) Q = 100 KJ / mol. 
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Fig. 16 The creep strains and stress relaxation in the ceramic coating as a function of time 
for the case of a lower heat flux 0.2 MW/m 2 . 
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Crack Initiation During Thermal Fatigue Tests 

The plasma sprayed ZTO 2 -Y 2 O 3 ceramic coatings contain microcrack networks 
with a typical crack width around 0.5-1 pm after processing. Therefore, initiation of larger 
cracks at the coating surface during thermal fatigue testing will not be a difficult process. 
The mechanisms of the crack initiation can be surface tensile stress induced cracking 
during cooling, and/or HCF peak compressive stress induced cracking at the heating stage. 
The surface tensile stresses are mainly generated by coating shrinkage after cooling due to 
the coating sintering and creep at temperatures. The pulsed laser induced temperature swing 
can generate locally high compressive stresses that could result in the surface coating 
fracture in a short time period. Since the laser HCF component will promote both the 
coating surface creep and the coating surface compressive cracking, the accelerated crack 
initiation and higher surface crack density at the coating surfaces are expected. This has 
been confirmed by this experiment. 

Fatigue Behavior of Thick Thermal Barrier Coatings under Thermal Cyclic Loading 

The fatigue crack propagation rates in a ceramic material under cyclic loads can be 
written as [17 ‘ 19] 


^ = CCa* (Kma* ~ *min V = C^AK” (4) 

dN 

where C, m and p are material dependent constants, /f max and K min are the maximum and 

minimum stress intensity factors, and A K the stress intensity amplitude, of the crack. 
Under the condition that K mm equals zero, Equation (4) can be reduced to the conventional 

Paris law relationship ™ 


— = CAK* (5) 

dN 

where q = m + p. During a superimposed thermal LCF and HCF testing, the surface 
vertical crack growth can be generally induced by both LCF and HCF components, as 
illustrated in Figure 17. The crack growth rate with respect to LCF cycle number can thus 
be expressed as 



( 6 ) 
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where Q and C 2 are constants, N* HCF is the characteristic HCF cycle number, A K LCF and 
AK hcf ^ stress intensity factors of the crack under low cycle and high cycle loads, 
respectively. The stress intensity factors are functions of crack geometry, crack length and 
stress magnitudes. It can be seen that the crack propagation rate depends not only on 
coating properties, but also on LCF and HCF parameters which define the stress states and 
fatigue mechanisms. 


A 



Cycle number 


Fig. 17 Schematic diagram showing crack growth resulting from thermal LCF and HCF 
loads. 


Low Cycle Fatigue Mechanism 

Under the present test conditions, the oxidation of the bond coat and the substrate is 
not important because of the low interfacial temperatures and short testing times. 
Therefore, the low cycle fatigue mechanism is primarily associated with coating sintering 
and creep at high temperatures. The time and elastic stress dependent, non-elastic strains in 
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the ceramic coating will lead to a tensile stress state during cooling, as schematically shown 
in Fig 18. This LCF stress under biaxial condition can be written as 

h £ 

C LCF = | £p ( a th » ' ~7 ~~ ” ( 7 ) 

0 1 v c 

where e p (a th ,T,t) is the strain rate resulting from ceramic sintering or creep, as has been 
described by Equation (2). The bond coat and metal substrate creep is not considered 
because of the low temperatures at the interfaces during the thermal fatigue testing. The 
LCF stresses as a function of time and coating layer depth are illustrated in Figure 19. The 
mode I stress intensity amplitude for LCF crack growth can be written as 

LCF = Z * [&LCF ~ & th ] ’ V fl&CO (8) 

where Z is a geometry factor associated with the crack configuration. Assuming that the 
crack does not grow under the compressive thermal stress o th , the stress intensity will 
depend primarily on <J LCF and the crack length a(i) . Therefore, the LCF crack growth rate 
will increase with time because of the increased stress G LC f level and the crack length. 
However, due to the stress <J LC f distribution profiles in the coating and its interactions 
with the ceramic/bond coat interface, the crack growth rate becomes more difficult to 
predict when the crack approaches the interface. Further work is underway to improve the 
understanding of the crack propagation and interface delamination. From Equation (8), it 
can also be expected that a faster crack growth rate will result with faster coating sintering 
and creep rates in the coating. 
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Fig. 18 


Sintering and creep at high temperatures 


tfLCF 



<3lcf 


Ceramic sintering and creep result in non-elastic strains (shown in shadowed 
area) at temperature, thus generating tensile stresses upon cooling. 


Time, hours 



Distance from the surface, mm 

Fig. 19 Tensile stresses are generated in the ceramic coating during cooling as a function 
of time and coating layer depth. These stresses are considered as a primary 
mechanism for LCF crack growth. 
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High Cycle Fatigue Mechanism 

The high cycle fatigue is associated with the cyclic stresses originated from the high 
frequency temperature fluctuation at the ceramic coating surface. Because this temperature 
swing results in significant thermal strains, considerable stresses will develop at the coating 
surface. HCF stresses are dynamic in nature, with a very short interaction time; therefore, 
stress relaxation can be neglected. The HCF stress amplitude is dependent on the 
temperature swing, and a stress level of 100 MPa can be induced at the surface by a 
temperature change of 250°C. With a surface crack in the coating, the HCF thermal loads 
can be equivalently acting on the crack by a wedging process, as schematically illustrated in 
Figure 20 (a) and (b). This wedging process, which provides an intrinsic mechanism for 
the HCF phenomenon, can be further enhanced by crack face shifting and spalled particle 
intruding, as shown in Figure 20 (c) and (d). Since the minimum HCF stress intensity 
factor equals zero, the net mode I stress intensity amplitude for this case can be expressed 
as ™ 


2 • P 1 + f(i) j TTT /n . 

AKihcf ~ / ' 9 V ^ ‘ a ( l ) (^ a ) 

and 

P - &HCF ’ (^b) 

where P is a concentrated load per unit thickness acting on the crack, b t is the load acting 

distance from the surface which is taken as laser interaction depth in the present study, 
o HCF is the HCF stress, a(i) is the crack length at the ith cycle. f(i) is a geometry factor, 
which can be related to the crack length a(i) and the interaction depth b i in the following 
form [21) 


/(0 = 


(A.Y 


n r 


0.2945-0.3912 


a(i) 


+ 0.7685- 


(AY 


-0.9942 


( b, f ( b ^ 


U(oJ 


+ 0.5094 




( 10 ) 


Note from the above that the stress intensity increases, in a linear manner, with 
increasing HCF stress o HCF and, by a more complicated function, with increasing 
interaction depth b t . The HCF stress is affected by the temperature swing AT, the thermal 
expansion coefficient a c and the elastic modulus E c of the ceramic coating. Figures 21-24 
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illustrate the relationship between the stress intensity factor and the normalized crack 
length, with various values of b t , AT, a c and E c of the coating. The results show that the 

stress intensity factor, thus the high cycle fatigue effect, decreases with increasing crack 
length, but increases with increasing the interaction depth, the temperature swing, the 
thermal expansion coefficient and Young's modulus of the ceramic coating. It should be 
noted that, depending on the coating stress state at high temperature, the HCF may affect 
crack propagation far beyond the laser interaction depth. This has been demonstrated in 
pure HCF cycling where high temperature swings, and therefore high thermal loads, were 
generated near the surface of the ceramic coating while the interior of the specimen 
remained cool ^ 2 \ This test condition was shown to cause not only surface crack initiation 
but also propagation deep into the coating, as shown in Figure 25. In fact, some of the 
cracks have reached the ceramic/bond coat interface after 5000 cycles when surface 
temperature swing was 700°C. In another experiment with a lower temperature swing 
from lower laser energy input, the crack growth was slower. 
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ceramic 


substrate' 



(b) crack growth under HCF conditions by surface wedging mechanism 




(c) enhanced surface wedging damage (d) enhanced surface wedging damage 
by surface crack face shifting by spalled particle intruding 

Fig. 20 Schematic diagram illustrating surface wedging mechanism during high cycle 

fatigue process. 
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Fig. 21 The relationship between the stress intensity factor amplitude and the laser 
interacting depth as a function of the normalized crack length. 
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Fig. 22 The relationship between the stress intensity factor amplitude and the temperature 
swing as a function of the normalized crack length. 
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Fig. 23 The relationship between the stress intensity factor amplitude and the thermal 
expansion coefficient of the coating as a function of the normalized crack length. 
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Fig. 24 The relationship between the stress intensity factor amplitude and the elastic 
modulus of the coating as a function of the normalized crack length. 
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Fig. 25 


A surface crack propagated deeply into the ceramic coating after 5000 
thermal shock cycles at a temperature swing of 700 °C. Each laser pulse 
heating and cooling cycle consisted of 0.1 second heating and 60 second 
cooling, respectively with interaction depth about 0.3 mm. 
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The Interactions between LCF and HCF Crack Growth 

Strong interactions between LCF and HCF have been confirmed by the prehminary 
experiments. More severe coating damage has been observed for the combined LCF and 
HCF tests compared to the pure LCF test. The higher crack density and wider crack width 
observed in the specimens with the higher test temperature suggest that significant coating 
sintering and creep are detrimental to the coating fatigue resistance. Higher heat flux near 
the beam center, as imposed by the spatial energy distribution of the Gaussian laser beam, 
resulted in increased surface cracking and spallation. This result can be expected because a 
higher heat flux will lead to not only a higher surface temperature and temperature gradient 
across the coating (generating more significant stress relaxation at temperature and thus 
more severe LCF damage after cooling), but also a greater temperature swing that enhances 
HCF failure. It seems to be true that both LCF and HCF are affected by the coating system 
configurations. In one dimensional coating systems such as angle iron comers and 
specimen edges, the cracking is less likely to occur in the less constrained direction, which 
is perpendicular to the one dimensional line direction. This result can be explained by 
considering that both the LCF and HCF stresses would be much lower in the less 
constrained direction. This implies that if a perfect bond coat strain isolation can be 
achieved, the coating fatigue resistance could be greatly improved. Further studies are 
required to obtain a better understanding of this phenomenon. 

The interactions between LCF and HCF lead to an earlier failure of the coating. The 
high cycle fatigue component promotes surface crack initiation and increases surface crack 
densities. This causes fast initial crack propagation near the coating surface according to the 
surface wedging mechanism, because of the extremely high stress intensity values at the 
initial stage. The longer cracks then increase the subsequent LCF stress intensity 
amplitudes, thus leading to a faster crack growth rate under the LCF mechanism. The LCF 
component will accelerate the subsequent HCF crack growth at high temperatures by 
predominantly two mechanisms. First, stress relaxation at high temperatures, which results 
from coating sintering and creep under LCF cycling, as well as from LCF induced crack 
formation and propagation, can significantly reduce the effective compressive stresses in 
the coating. The HCF crack growth will be facilitated by this process. Second, the coating 
surface sintering under LCF cycles could considerably increase the coating elastic modulus. 
A higher coating modulus will lead to higher HCF stresses, and thus enhance the HCF 
crack growth. 
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SUMMARY 


The present analysis presents a detailed explanation of the processes occurring 
when thick TBCs are subjected to combined thermal low cycle and high cycle fatigue. This 
work also provides a framework for developing strategies to manage ceramic layer 
sintering and creep, thermal expansion mismatch, and other characteristics of the coating 
system. For example, since ceramic sintering and creep are detrimental to the coating 
system, creep resistant coatings, especially near the surface region, would be expected to 
improve the durability of the system. In addition, since it is well known that LCF behavior 
is closely related to the thermal expansion mismatch, good strain isolation provided by 
well-structured and compliant bond coats would further improve the fatigue resistance. The 
HCF resistance could be effectively achieved by high compressive stresses in the coating. 
A high toughness top surface layer with low modulus and thermal expansion coefficient 
would also improve the HCF fatigue life. The relative importance of LCF and HCF crack 
growth and its correlation with coating stress states, creep behavior and bond coat 
properties need to be carefully investigated. 

CONCLUSIONS 

1. Strong interactions between LCF and HCF have been observed in preliminary 
experiments. The combined LCF and HCF tests induced more severe coating damage 
compared to the pure LCF test. Significant coating sintering and creep, which are 
confirmed to accelerate both the LCF and HCF failure by experiments, are detrimental to 
the coating fatigue resistance. 

2. In the absence of severe interfacial oxidation, the LCF mechanism is closely related 
to coating sintering and creep phenomena at high temperatures. The stress relaxation, at 
temperature, induces tensile stresses in the coating after cooling. However, the HCF 
mechanism is associated with the surface wedging process. The HCF damaging effect will 
increase with the heat flux, thus the temperature swing, the thermal expansion coefficient 
and the elastic modulus of the ceramic coating, as well as the HCF interaction depth. 
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APPENDIX NOMENCLATURE 


t c , t b and t s 
a 

I 0 and I (r) 
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o , <J and 

^ jtotal 

(J th and 
r and AT 
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t and tj 
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Q 

e p 

K K 

Jv max ’ Jv min 

and A K 


Ceramic coating, bond coat and substrate thicknesses, mm 
Thermal expansion coefficient, m/m-°K 

Laser irradiance or power density at the beam center and distance r 
from the center , MW/m 2 

Laser beam total power, W 

Thermal conductivity, W/m-°K 

Density, kg/m 3 

Heat capacity, J/kg-K 

Young's Modulus, GPa 

Poisson's ratio 

Thermal stresses, residual stresses and total stresses in coating 
systems, MPa 

Thermal stress and initial thermal stress in ceramic coating, MPa 

Temperature and temperature swing, °K 
Gas constant, J/mol-°K 
Time, sec. 

Pre-exponential constant for ceramic coating creep 
Stress exponent for ceramic coating creep 
Time exponent for ceramic coating creep 
Activation energy for ceramic coating creep, J/mol 
Ceramic coating creep strain rate, 1/sec 
Ceramic coating creep strain 

Maximum and minimum stress intensity factors, and the stress 

1/2 

intensity amplitude, of the crack, MPa • m 
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LCF an ^ 

^K\hcf 
C, Cj and C 2 
ra, /? and q 

N and N* HCF 
a(i ) 
h 

P 

Z and f(i) 


Mode I stress intensity factor amplitudes of the crack under low cycle 

i /2 

and high cycle loads respectively, MPa • m 
Constants 

Stress intensity exponents in fatigue, and q = m + p 
LCF cycle number and HCF characteristic cycle number 

Crack length at the ith cycle, mm 
Laser interaction depth, mm 

Concentrated load per unit thickness acting on the crack, N/m, 
P = °HCF ‘ 

Coefficients associated with the crack configuration 
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expansion coefficient and the elastic modulus of the ceramic coating, as well as the HCF interacting depth. A good 
agreement has been found between the analysis and experimental evidence. 
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